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ABSTRACT 

Context. The identification of long-gamma-ray-bursts (LGRBs) is still uncertain, although the coUapsar engine of fast- 
rotating massive stars is gaining a strong consensus. 

Aims. We propose that low-metallicity Be and Oe stars, which are massive fast rotators, as potential LGRBs progenitors. 
Methods. We checked this hypothesis by 1) testing the global specific angular momentum of Oe/Be stars in the ZAMS 
with the SMC metallicity, 2) comparing the ZAMS {Q./Q.C, M /Mq) parameters of these stars with the area predicted 
theoretically for progenitors with metallicity Z = 0.002, and 3) calculating the expected rate of LGRBs/year/galaxy and 
comparing them with the observed ones. To this end, we determined the ZAMS linear and angular rotational velocities 
for SMC Be and Oe stars using the observed V^sini parameters, corrected from the underestimation induced by the 
gravitational darkening effect. 

Results. The angular velocities of SMC Oe/Be stars are on average (Sl/fJc) =0.95 in the ZAMS. These velocities are 
in the area theoretically predicted for the LGRBs progenitors. We estimated the yearly rate per galaxy of LGRBs and 
the number of LGRBs produced in the local Universe up to z—0.2. We have considered that the mass range of LGRB 
progenitors corresponds to stars hotter than spectral types BO-Bl and used individual beaming angles from 5 to 15°. 
We thus obtain -R£q^b ~ 10~^ to ~ 10~^ LGRBs/year/galaxy, which represents on average 2 to 14 LGRB predicted 
events in the local Universe during the past 11 years. The predicted rates could widely surpass the observed ones [(0.2- 
3)xl0~^ LGRBs/year/galaxy; 8 LGRBs observed in the local Universe during the last 11 years] if the stellar counts 
were made from the spectral type B1-B2, in accordance with the expected apparent spectral types of the appropriate 
massive fast rotators. 

Conclusions. We conclude that the massive Be/Oe stars with SMC metallicity could be LGRBs progenitors. Nevertheless, 
other SMC O/B stars without emission lines, which have high enough specific angular momentum, can enhance the 
predicted -Rlgrb rate. 

Key words. Gamma rays: bursts - Stars: early-type - Stars: emission-line. Be - Stars: fundamental parameters - 
Galaxies: Magellanic Clouds 

1. Introduction star scenario (jTutukov &: Fedorovall2007t ). However, proper 

identification of both the progenitor and the final nature of 
Since their discovery, data on gamma ray bursts (GRBs) the phenomenon are still fairly uncertain. Nevertheless, a 
considerably increased, which today allow us to aim at pre- kj^d of consensus is gaining in the community that the col- 
; liminary statistical conclusions. The properties of GRBs lapsar en gine during a su pernova SNIb,c explosion, as pro- 
and their possible connection to SN were reviewed by posed bv lWooslevI (fl99l . might lead to allow the explana- 
IWoosley fc Bloom| ([20061) and [Fryer et al.| ([20071). Thus, tion of LGRBs. According to this hypothesis, the explosion 
there are at least 2 classes of GRBs according to the dura- follows a massive star collapse to a black hole. The circum- 
tion of the phenomenon: 1) short bursts that last less than stellar dis k accretes onto the black hole and a bi-polar jet 
1 s and 2) long bursts, which are longer than 1-3 s here- is formed (iHirschi et al.l[2oH5h . The infalling material must 
after LGRBs. Although the short GRBs might correspond have enough angular momentum to remain in a disc before 
to a violent merging of two compact objects, evidence is accretion. 

growing that suggests that the disrupti ongf a massive star observational findings and the latest theoretical 

can be behmd an LGRB jPryer et al. |M- A possible, developments are steadily giving best founded constraints 
but rare, thn-d class of GRBs can be explamed by a bmary- understand t he LGRBs phen omenon. The observations 

carried out by llwamoto et al.l (|l998l l2000t) support the 

Send offprint requests to: C. Martayan idea that mass ive fast-rotat i ng st ars are at the origin 

Correspondence to: Christophe.Martayan@eso.org of the LGRBs. iThone et all ([2008D find that the LGRB 
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GRB060505 is hosted in a low-nietallicity galaxy, which 
is characterized by a high star-formation rate. The event 
seems to come from a young environment (6 Myrs) and from 
an object about 32 Mq . The latest models have tried to re- 
produce the evolution of stars until the end of their lives. 
In particular, for massive stars, they have been worked out 
until the GRBs phase. In the search for SNIb,c progenitors, 
WR stars with He-rich envelopes are recognized as pos sibly 
behind the GRBs events. In fact, [Hammer eTall (pOOeh find 
that the LGRBs occur in areas of galaxies with WR stars. 

Rotation has been recognized as a key point for 
understanding the appearance of GRBs (iWooslevI [19931 
iHirschi et aI1l2005t [Voon et al.ll2006[) . Accordingly, to keep 
a large amount of angular momentum up to the last evo- 
lutionary phases before the collapse, GRBs progenitors 
should be massive objects with low i nitial metallicities 
and p ossibly displa y anisotropi c winds (jMevnet fc MaedeH 
120071) . Thus, from lYoon elall (|2006[ ) it seems that WR 
stars with metallicities Z < 0.002 can be p rogenitors of 
GRBs. In the same sense, iModiaz et all (120081) find that the 
SNIc-GRB association occurs when metallicities are lower 
than 0.6 x Zq. Becau se stars with l ow me tallic ities can on 
averag e rotate faster, iHirschi et al.l (l2005h and lYoon et aD 
(|2006[ ) foresee that the number of GRBs must grow as the 
redshift increases, which seems to be confirmed observa- 
tionally (jFrver et al.ll2007t) . It might stiU be that the first 
massive stars, which were very metal-poor stars, were fast 
rotators. This is an additional reason for the frequency of 
GRBs being higher as the redshift increases. In this sense, 
iKewlev et al.l (|2007t ) find a link between the location of cos- 
mological LGRBs and the very low-metallicity galaxies. 

Thanks to fast rotation and the concomitant eSicient 
mixing of chemical elements, massive stars can undergo 
quasi-chemically homogeneous evolution to end up as he- 
lium WR stars satisfying the requirements for the co l- 
lapsar scenario (lYoon et all l2006t Ivan Marie et al.ll2008D . 
lYoon et all (|2006[ ) have calculated the quasi-chemically ho- 



mogeneous evolution of magnetized massive stars. They 
produced diagrams of LGRBs progenitors as a function of 
their ZAMS rotational velocities and masses and of different 
initial metallicities. From these diagrams, it emerges that 
there must be an upper limit to the initial metallicity of 
LGRBs progenitors, which approximatively corresponds to 
that of the Small M agellanic Clouds (SMC). According to 
iMaeder et al.l (|1999L and references therein), t he SMC av- 
erage metallicity is Z~ 0.002. The diagram of (jYoon et al.l 
|2006[) also indicates that, for metallicities Z < 0.002, the 
WR phenomenon can appear in stars having lower masses 
than those in the Milky W ay (hereafter MW). In this con- 
text, iMartins et all (|2009D have observed several SMC WR 
stars, whose evolutionary status and chemical properties 
can be understood if they are fast rotators. 

According to the listed observational and theoretically 
inferred requirements, stars might be potential LGRBs pro- 
genitors if they 

— are massive enough at the end of their evolution to form 
a black hole; 

— have lost their hydrogen envelope and have a fast- 
rotating core; 

— are formed in low-metallicity environments, where there 
must also be high star-formation rates to ensure having 
enough massive stars. 



Since the lower the metallicity the faster the rotation 
and the lower the mas s of stars that c an undergo the 
WR phase, according to lYoon et al.l (|2006l) . the most mas- 
sive B-type stars, as well as 0-type stars, could become 
WR stars if they rotated fast enough. The required ro- 
tational velocities must be close to those what are typ- 
ical of Oe/Be stars. Then, we simply ask whether the 
more massive stars displaying the Be phenomenon today 
in the SMC, can be potential progenitor s of LGRBs. This 
possib ility has already been raised by IWooslev fc Hegerl 
(|2006t ). We recall that the "Be phenomenon" appears in 
the main sequence evolutionary phase of fast-rotating O- 
and B-type stars. They are characterized by emission lines 
produced in a decretion disk formed by continuous and 
episodic matter ejections from the central star. These ob- 
jects are the fastest known rotators in the main sequence, 
and their rotation can b e closer to the critical one when 
the metallicity is lower (iDomiciano de Souza et "ahl 2003[: 
Meilland erall 120071: iFremat et al.l 120051: IVinicius et"al] 



20061: iMartavan et al.ll20oW 



Our research on whether the SMC Oc/Be stars can be 
progenitors of LGRBs is based on three tests: 1) deter- 
mination of the global specific stellar angular momentum 
(Sect. [5]), 2) inference of the ZAMS angular velocity ra- 
tios of the studied SMC Oe/Be stars and comparison with 
the model predicted requirements to be LGRB-progenitor 
(Sect. [3]), and 3) estimation of the LGRBs rate based on the 
SMC Oe/Be- and fast-rotator population, and its compar- 
ison with the observed ones (Sect.H]). A general discussion 
is given in Sect. |5j 



2. Inference of the ZAMS rotational velocities 

2.1. Main sequence linear equatorial rotational velocities 

The results presented in this section are based on the 
FLAMES-G IRAFFE (iPasauini et aD I2002D observations 
discussed in lMartavan et al.l (|2007f I]~To the previous sain - 
ple of 131 Oe/Be stars analyzed bv lMartavan et all (|2007t ). 
we added here the most massive Be stars and Oe stars found 
in the SMC. The projected rotational velocity Vsin i of each 
star w as determined with the GIRFIT code (jFremat et al.l 
120061) . The average values of Vsin i by mass category thus 



obtained are given in Table [T] We know, however, that 
Oe/Be stars are fast rotators, and the gravitational dark- 
ening eff'ect (hereafter CD) induces systematic underesti- 
mations of the Vs in i parameters (iTownsend et al.l 12004 
iFremat et al.l l2005t ) . In this work we have corrected the 
measure d Vsini quantities using similar curves to those 
given m IFremat et all (|2005D but recalculated for the SMC 
metallicity. The corrected average rotational parameters 
per mass category are presented in Table [T] as (ysini)corr- 
Figure [T] shows the (ysini)corr per mass category of the 
studied SMC Oe/Be stars, compared with three new theo- 
retical tr acks of true eq u atorial rotational velocities calcu- 
lated bv lEkstrom et al.l (l2008bD for metallicity Z= 0.002, 
which suits the SMC cluster NGC330 and its environment. 
In what follows, all models for SMC stars are for Z= 0.002. 
To compare the theory with observations, all model ve- 
locities were multiplied by 7r/4 = (sini) to simulate the 



^ Based on observations at the European Southern 
Observatory, Chile under project number 072.D-0245(A) 
and (C). 
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rand om distribution of the inclinati on i of the rotational 
axis (jChandrasekhar fc MiinchlllQSClf ). Once the fundamen- 
tal parameters of all studied stars had been determined, 
the models of rotating stars were chosen so as to match the 
average true rotational velocity (i.e. V = (4/7r)ysini) per 
mass category at their respective average age. The average 
true rotational velocities, V, per mass category have uncer- 
tainties ranging from 20 to some 40 km s~-'^(cf. Table [T]). 
The inferred average ZAMS velocities for these mass cat- 
egories are then affected by uncertainties that are of the 
sam e order of magnitude, i.e. SV ~ ±30 km s^^, on aver- 
age (jMartavan et al.l2007l their Fig. 9). For the near critical 
linear velocities, V ^ Vc, and stellar masses M > 17 Mq, 
the corresponding jiJl^/Vcl translates into |5(ri/ric)| ^ 0.05 
as the average uncertainty of the derived angle velocity ra- 
tios {Vc and fic are the linear and angular critical velocities, 
respectively) . 

The obtained today Jl/rJc ratios for all mass groups arc 
in the interval 0.99 S^/^cS l-O- Fig- H shows the loci of 
today (Vsmi) against the mass corresponding to il/ilc = 
0.99 and 1.0. For comparison, this figure also includes the 
(Vsini) curve for n/n^^OM. 



A similar comparison in iMartavan et al.l ([2003) , made 
with observed rotational velocities that were not cor- 
rected for the underestimations induced by the GD ef- 
fect, suggestes that the SMC Be stars rotate on average 
at fi/f2c=0.95. Within the observational uncertainties, the 
corrected velocities for the GD effect shown in Fig. [T] likely 
correspond to n/rtc~0.99. We note that r2/f2c=0.99 is for 
a linear velocity ratio V/Vc = 0.97 or for a ratio of surface 
centrifugal to gravitational accelerations 77 = 0.90. 



Table 1. Values of the rotational velocities for Be and Oc 
stars for several average mass categories 



Mass category 


A 


B 


C 


D 


E 


{M/Mq) 


3.7 


7.6 


10.8 


13.5 


23.7 



N 


14 


82 


13 


15 


7 


(Vsini) 


277 


297 


335 


336 


420 


(ysini)corr 


300 


326 


376 


377 


457 


log[age(yr)] 


8.3 


7.6 


7.3 


7.2 


6.9 




400 


437 


513 


514 




/-i/\ncw 
V /ZAMS 


448 


487 


562 


563 


683 


^VzAMS 


35 


27 


23 


23 


44 


/T/X magnet 
V /ZAMS 


389 


423 


488 


489 


643 



N = number of stars by mass category 

(^)zAMS ~ average true rotational velocities in km/s 

without correction for GD fMartavan et al. 2007) 

(^)zAMS = average true rotational velocities in km/s 

of the enlarged stellar sample, corrected for GD effect 

evz^^fs^= error of (^)zAMS in km/s 

(^)zAMS = average true rotational velocities in km/s 

of the enlarged stellar sample, 

corrected for models including magnetic fields 



2.2. ZAMS rotational velocities 



Using the same approach as in 

IMartavan et al.l (I2007D. whi c h was base d on models calcu- 
lated bv iMaeder fc Mevnetl (|2001f ) and iMevnet fc Maedeil 
(120021) . we re-determined the ZAMS rotational velocities 
of the SMC Oe/Be stars with a new and more complete 
grid of curves representin g the evolution of rota tional ve- 
locities at low metallicity ( Ekstrom et al.l[2008bf ). The ob- 
tained results are given in Table [1] For the sake of com- 
parison, this table also gives the previous average values of 
the ZAMS rotational velocities. The new results are shown 
in Fig. [2j where t he theoretical Z AMS r otational veloci- 
ties calculated by lEkstrom et al.l ()2008b() are also drawn 
for the SMC metallicity and for the angular velocity ra- 
tios 17/17c=0.90, rj/rJc=0.99 and ^1/^^=1.00. We notice 
then that SMC Oe/Bc stars in the ZAMS rotate on aver- 
age with fi/f2c— 0.95. The models do not include magnetic 
fields, thus the mean rotational velocity on the MS is lower 
than it would be for magnetic models. We might thus over- 
estimate the initial ratio ft/ftc, which should be considered 
here as an upper limit. To take the difference in the models 
into account with magnetic fields, we compared the val- 
ues of ZAMS velocities for g i ven ^l/ ftr, values between the 
model s from lEkstrom et al.l ()2008b[ ) and from lYoon et al.l 
(|2006D . The values oiYzAMS in the last row of Table |T] are 
then corrected by the corresponding scale factor. 
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Fig. 1. Comparison of the average observed rotational ve- 
locities SMC of Oc/Bc stars with the locus of theoreti- 
cal velocities for the SMC metallicity. The curves are for 
model stars that attain the surface angular velocity ratios 
f^/^2c=0.90, ri/^2c=0.99 and fi/fJc=l-00, in the MS evo- 
lutionary phase. 



The trend of points inferred with the observation in 
Fig. [2] can be represented with (V^|^ms>- 301xM°-26±° °5 
km s~^. For the inference of the specific angular mo- 
mentum that must have the circumstellar disks hanging 
over the collapsed stellar cores, it is instructive to esti- 
mate the average specific angular momentum, ( J/M) zams^ 
that the SMC Oe/Be stars have in the ZAMS. To this 
end, let us write {J/M)zams = kG{V§^]^jg){RzAMs), 
where fcc is the gyration radius and Rzams the stellar 
radius of stars in the ZAMS. The fee and Rzams radii 
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suited to ri/Oc = 0.95, chosen according to what is sug- 
gested by the observed points in Fig. [2j varies with the 
mass as kc = 0.017 x A/"-^'^^^ '-'^, while the radius varies 
as iR7.AMfi)= 1.72xM°-^^ ±o.02^ (IZoredll98fil: IZorec et al.l 
Il988t lEkstrom et"alll2008bf ). By taking the decrease in the 
stellar moment of inertia induced by a rigid rotation at 
fi/r^c = 0.95 into account, we obtain the following ex- 
pression for {J/M)zAMS, valid for the whole mass interval 
3 < M/Mq < 30; 



(J/M) 



ZAA/s - 6.0±0.1 X lO'^f— — ) cm^s \(1) 



which means that the average specific angular momentum 
of Oe/Be stars of the SMC in the ZAMS, considered as rigid 
rotators, increases linearly with the mass. Using moments 
of ine rtia calculated with stellar models at rest, iKawaleij 
(|l987t) obtained for a mixed sample of B-|-Be stars in the 
MW that (J/M) ~ Mi-43±o i6. Since the specific angular 
momentum of the last stable orbit aroun d a Sch warzschild 
black hole scales with the mass M/Mq as (|Yoon et al.„2006h 
(3 times more than for a Kerr black hole) 
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Fig. 2. Inferred ZAMS rotational velocities of SMC Oe/Be 
stars (triangles) compar ed with model SMC Z AMS ro- 
tational velocities from lEkstrom et all ()2008b[ ) (dashed 
curves) . 



(J/A/)5BH = 1.5x 10i«(^ 



? — 1 
cm s 



(2) 



the SMC Oe/Be stars would then have a specific angu- 
lar momentum only a few times greater than the one 
that is supposedly appropriate for producing a stable disk 
over the black hole after the stellar collapse. Since part 
of this angular momentum will be lost during the evolu- 
tionary following the MS phase, the SMC Oe/Be stars do 
not seem to have angular momenta that is to prevent the 
formation of circumstellar dis ks and/or the expected jets 
( MacFadven fc Wooslevlll999t) . 

Because it may affect the properties of Oe/Be stars, 
after we take the average of the fi/f^c rates over the masses 
in the ZAMS and in the MS, the equatorial velocity of these 
objects decreases from the ZAMS to the MS by a factor 
0.76, while n/Q.^ varies from 17/f7c = 0.95 to fi/f]c = 0.99. 
This change in fi/fic implies that the ratio V/Vc varies 
from 0.87 to 0.97 and the ratio 77 of the surface equatorial 
centrifugal to gravitational acceleration passes from 0.76 to 
0.90. 



3. LGRBs progenitors and the Oe/Be stars 

3.1. ZAMS rotational velocities of LGRBs progenitors 
compared to those of Oe/Be stars 

Thanks to the combined effects mainly of rotation and 
metallicity, a series of models suggest that massive fast- 
rotating stars can fulfill the requirements imposed by 
the coUapsar engine and be so pr ogenitors of LGRBs 
(jHirschi et alJl2005HYoon et al.ll2006D . To comply with the 
SN-LGRB relation, the involved stellar masses need to 
be higher than 10 Mg. Thus, while fllirschi et al. (2005i) 
conclude, on the basis of non magnetic models, that for 
the SMC metallicity stars with masses between 32 and 
92 Mp) are able t o produce LGRBs, the latest models by 
lYoon et al.l (|2006t ). which include magnetic fields decrease 
the lower limit of masses to those of early B-type stars. 
According to these models, fast rotation favors a strong 
mixing of chemicals, which means that the stars can evolve 



as quasi-chemically homogeneous and finally end up as mas- 
sive helium stars. Since the lower the metallicity, the faster 
the stellar rotation can be, these phenomenon can in princi- 
ple be favored in media with low metallicity, which accord- 
ingly impede strong stellar winds. In turn, this makes the 
stars cons erve high amounts of angular momentum. The 
models by lYoon et all (|2006t ) predict a lower limit to the 
ratio of equatorial linear velocities V/Vc (Vc is the critical or 
breakup velocity) as a function of mass for the LGRB pro- 
genitors that do not seem to change strongly according to 
the metallicity, but the limiting V/Vc ratio increases if the 
stellar mass is lower. Moreover, the zone above the thresh- 
old V/Vc ratio of potential LGRB progenitors widens to- 
wa rds lower and higher ni asses as the metallicity decreases. 

ICampana et al.l ()2008l and references therein) indicate 
that the LGRB GRB060218 had a progenitor with an ini- 
tial mass of 20 Mq, i.e., a BO-Bl star, where the metallicity 
is Z = 0.004, similar to the SMC. The chemical composition 
of the nebulae surrounding the star, if it is assumed that it 
reflects the chemical abundance of the progenitor, indicates 
that the progenitor was a near critical rotator, as in the "Be 
phenomenon". Such a fast rotation is required to produce 
an efficient mixing process in the star, which favors the stel- 
lar chemically homogeneous evolution and makes it possible 
for the environing nebulae reflect the stellar abundances. 

To compare the models with the inferred ZAMS ratios 
fl/ric, we transformed the 'observedQ (T^/Vi;)zAMS ratios 
of the Oe/Be star subsamples in the SMC into fl/ilc ratios 
by considering the la dispersion of masses and equatorial 
velocities intervening in each averag e. The V/Vc and f^/f^c 
ratios for massive stars are related as ([Chauville et al.l2001t) 



" 1.38lf [l-0.276(f)^] 



(3) 



which is also valid for the SMC metallicity. These results are 
given in Tabled] The indicated mass category, i.e. AztcrM, 



^ Actually, (V/V'c)zams must be understood as values inferred 
from observed V/Vc ratios through model predicted tracks of the 
stellar rotational velocity. 
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should be understood as A'/(A) ± ctM. The boldfaced value 
of the f2/r2c in % corresponds to the average mass (M) of 
the corresponding subsamplc. The natural upper limit of 
the surface stellar rotation is obviously given by V jV^ = 
1.0. In Table [5] all values V jV^ > 1.0 are the arithmetic 
consequences of the imposed symmetric la dispersions and 
they should be understood as V/Vc — >■ 1.0. In Tabled the 
values of Table [2] are corre c ted fo r the effect of magnetic 
fields from the lYoon et al.l (|2006D models as explained in 
Sect. O 

Table 2. Ratios of the ZAMS per mass category of 

SMC Oc/Be stars with the boldfaced values corresponding 
to the average Q/flc ratios of subsamples. 



Mass category 








A 


0.89 


0.79 


1.01 


A-ctM 


0.94 


0.83 


1.04 


A+ctM 


0.85 


0.75 


0.96 


B 


0.87 


0.68 


1.08 


B-o-M 


0.93 


0.73 


1.15 


B+o-M 


0.82 


0.64 


1.01 


C 


0.92 


0.71 


1.08 


C-ctM 


0.94 


0.72 


1.09 


C+ctM 


0.91 


0.70 


1.06 


D 


0.91 


0.63 


1.07 


D-ctM 


0.95 


0.66 


1.13 


D+crM 


0.87 


0.60 


1.02 


E 


0.92 


0.78 


1.02 


E-(tM 


1.03 


0.88 


1.15 


E+(tM 


0.82 


0.70 


0.91 



The inferred ZAMS zones of n/ftc ratios per mass cat- 
egory of SMC Oc/Bc stars given in Table |3] are compared 
in Fig. [3] with the theoretically predicte d region containin g 
the LGRB progenitors for Z = 0.002 bv lYoon et all (|2006D . 
One can see that only the more massive Oe/Be members, 
found in the "E-category" , may re ach the zone of pot ential 
LGRBs progenitors. According to lYoon et al.l (|2006[ ). the 
ratio of the core to envelope mass is low in the less-massive 
0/B stars, which favors an efficient braking of their core 
and thus a significant reduction of its specific angular mo- 
mentum. These stars were then not included in our counting 
as potential LGRBs progenitors. Nevertheless, among the 
less massive OB objects, there may be many that are pos- 
sibly storing high amounts of angular momentum, thanks 
to strong internal differential rotation. In these objects, the 
core specific angular momentum can then be higher than 
the one predicted with stellar evolutionary models con- 
structed with an initial upper limit of kinetic energy im- 
posed by the rigid rotation in the ZAMS. Thus, in spite of 
the limiting masses imposed by the above mentioned the- 
oretical predictions, many more stars could be potential 
LGRBs progenitors. The same arguments can also be ex- 
tended to stars, of the same and other masses, which are 
here excluded from the counting because of their apparent 
low surface rotation. We discuss this issue in Sect. [5j 

In Fig. 13] the shaded "E" zone seems to show somewhat 
an unusual behavior compared to other mass categories. 



Table 3. Ratios of the ZAMS n/ftc per mass category of 
SMC Oe/Be stars corrected for the effect of magnetic field 
from the compari son between Ekstrom et al. (2008b) and 
lYoon et al.l (l2006l ) models (boldfaced values correspond to 
the average il/flc ratios of subsamples). 



Mass category 








A 


0.77 


0.69 


0.88 


A— <tM 


0.82 


0.72 


0.90 


A+aM 


0.74 


0.65 


0.83 


B 


0.76 


0.59 


0.94 


B— ctM 


0.81 


0.63 


0.99 


B-^ctM 


0.71 


0.56 


0.88 


C 


0.80 


0.62 


0.94 


C-crM 


0.82 


0.62 


0.95 


C+aM 


0.79 


0.61 


0.92 


D 


0.79 


0.55 


0.93 


D-ctM 


0.82 


0.57 


0.98 


D+(tM 


0.76 


0.52 


0.89 


E 


0.83 


0.70 


0.92 


E-crM 


0.93 


0.79 


1.03 


E+aM 


0.75 


0.63 


0.82 



Although we are not able to specify which of them is dom- 
inant, four reasons may be used together to explain it: 1) 
the smallness of the stellar sample which by chance perhaps 
favors the highest Vs'mi values; 2) fast evolution of mas- 
sive stars in low-metallicity regions that enables the initial 
high angular momentum to be retained; 3) the overesti- 
mated Vsmi of massive Oe/Be stars because of some sys- 
tematic non ro tational line broadening effect, i.e. elect ron 
scattering, etc. (jBallereau et al.lll995l:IZorec et al.lll992D : 4) 
the stronger mass-loss in the more massive stars, so that 
their surface D,/D,c ratio is le ss likely to increase durin g the 
MS lifes pan, see for ex amples 'Mevn et fc Maedeil (|2000L Fig. 
11) andlMaeder fc Mevnet (,200L Fig. 4). 



iHunter et al.l 



indicated that, according to the 



Vsini distributions of B stars they observed in NGC 
346, the areas of LGRBs progenitors c ould be more ex- 
tended than expected from lYoon et al.l (l2006l ) model pre- 
dictions. T he ag e of this cluster is roughly 3-5 Myrs 
(|Nota et all 120061: iBouret et al.ll2003l ). there have been re- 
cent episod es of star format i on an d it hosted at least 1 
supernova (|Gouliermis et al.l l2008[) . In this region there 
are also several classical Be stars and other emission-line 
objects (IWisniewski &: Biorkmaii 20061: IWisniewski et al.l 



120071: iHunter et al.ll2008t iMartavan et al.(l2010D T 



3.2. Evolution of Oe/Be and WR stars 

Some stars could start in the ZAMS as born Oe/Be objects, 
but lose this property during their MS life-span due to the 
spin-down produced by efficient mass- / angular momentum- 
loss phenomena. This can be the case of the massive Be 
stars in the MW, which a re only found in the first half of 
the MS (|Zorec et al.ll2005D . Contrary to the MW, the mas- 
sive Oe/Be stars and Oe stars in the SMC are also present 
in the second half of the MS. These stars can then continue 
their post-MS steps by keeping fast rotational rates and 
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Fig. 3. Comparison of the (0/Oc, M/Mq) parameters in 
the ZAMS per mass category for SMC Oe/Be stars with 
the area of long gamm a ray burst s prog enitors at SMC 
metallicity predicted bv lYoon et al.l (|2006[ ). 



evolve as chemically quasi- homogeneous ob jects towards 
the fast-rotating WR phase (lYoon et al.ll2006D . This p redic- 
tion seems to be supported bv iMartins et all (|2009h . who 
find that several SMC WR stars are certainly fast rotating 
objects. 

The finding in Sect. l3.11 that the more massive members 
of the SMC Oe/Be po pulation overlaps the LGRBs pre- 
dicted progenitor area (jYoon et al.ll2006l) . thanks to their 
masses and (ri/ric)zAMS rotation rates, which are today in 
the second half of their MS evolutionary phase, makes them 
designed candidates for such chemically quasi-homogeneous 
objects that can evolve towards the WR phase and become 
LGRB progenitors later on. 

iMartavan et al.l (120101 ) find that the "Be phenomenon" 
can be shared by more numerous stellar populations, fa- 
vored possibly by the lower SMC metallicity, than in the 
MW. Since the theoretical predictions show that the lower 
the metallicity the lower is the stellar mass that can enter 
the LGRB progenitor area, it is expected that among the 
massive stars of the first generations, which are very metal 
poor, there might be a significantly large number of objects 
that could have undergone the "Be phenomenon". These 
stars could then have had the required physical conditions 
for being potential LGRB progenitors. 



4. Predicted LGRBs rates from the Oe/Be star 
population and comparison with the observed 
ones 

Another test of the ability of the Oe/Be stars in the SMC 
to be LGRBs progenitors is based on the estimation of the 
expected rate of LGRB events that these stars can afford. 
In the less favorable case, the predicted frequency must be 
of the same order of magnitude as the observed ones for 
this stellar population be considered a plausible candidate. 



4.1. The SMC Oe/Be stellar population and selection of 
candidates 

To estimate the LGRB rate from the SMC Oe/Be steUar 
population, we need to identify the right spectral types to 
fill the mass requirements of stars entering the theoretically 
predicted LGRB progenitor zone shown in Fig. [31 To this 
end, and knowing that we are dealing with fast rotators, 
we used the models of stellar evolutio n with rotation for 
stars w ith SMC metallicity calculated bv lMaeder fc Mevnetl 
(|200lh . We interpolated the evolution paths from the ZAMS 
to the TAMS for the masses M = 18Mq and M = 33Mq. 
The (logL/L0,Teff) corners thus obtained in the theoreti- 
cal HR diagram were transformed into absolute magnitude- 
color pairs [My,{V — I)o] using standard bolometric cor- 
rections and the OGLE-III color-magnitude calibrations in 
current use. It follows from this that the main sequence 
visual absolute magnitudes of SMC stars with M = 18Mq 
range from My = — 2.8 mag to My = — 4.4, while those of 
33Mq do it from My = —3.7 mag to My = —6.3. A mass- 
spectral type corresponde nce for stars in the ZAM S has also 
been recently obtained bv lHuang fc GiesI (|2Q06( ). 

In accordance with the MK spectral-types for dwarf 
stars in the SMC calibr ation against the visu al absolute 
magnitude My used in (jMartavan et al.ll2Q10r) . the stars 
entering the LGRB progenitor zone in Fig. |3] correspond 
then to those labeled with B1-08 spectral types for rest 
and em ission-less stars. However, according to lHunter et al.l 
(|2008[ ) the mass of the LGRBs progenitors could be as low 
as 14Mo. In this case, one should count stars a little cooler 
than spectral type Bl, i.e. from B1.5, up to 08 as potential 
LGRBs progenitors with SMC metallicity. Stars from 14 to 
33 Mq evolve in the MS phase, so that their respective ab- 
solute magnitudes My remain rather constant and that the 
{V — /)o color changes by no more than 6{V — I)o ^ 0.04 
mag for 14 Mq to S{V - /)o 0.07 mag for 33 Mq. 
Moreover, in our [My, {V — /)o] diagram, all counted OB 
stars are located in quite a narrow strip of My a. s a func- 
tion of (y — /)o (see Fig. 5 in lMartavan et al.]|2010() . so that 
there cannot be much confusion regarding the identification 
of stars by their masses with the employed photometric cri- 
terion. 

The absolute magnitude intervals for 18 and 33M0 
stars were established from evolutionary models where the 
{log L/Lq, Toff) parameters are averaged over the rotation- 
ally deformed stellar surfaces. Actual stars show, however 
apparent hemisphere-dependent spectroscopic and photo- 
metric characteristics. There are then at least two reasons, 
rotation- and Be phenomenon-related, for which still cooler 
spectral types than Bl should enter the counting: 1) fast- 
rotating stars appear with cooler app a.rent spectral types 
than stars with the same mass at rest (|Fremat et al.ll2005l: 
IMartavan et al.l[2007f ): 2) the largest number of candidate 
stars are selected using photometric colors in the visible 
and near-IR, where both rotation and the presence of cir- 
cumstellar disks make the stars appear strongly reddened. 
We recall that reddening due to the circumstellar disks 
A{B — V) > 0.05 mag are currently observed in Oe/B e 
stars (jMouitahid et al.lll998l Il999t IMartavan et"al]l2010f ). 
These two effects make that the stars with photometrically 
determined spectral types can appear from two to three 
subspectral types cooler than expected for their intrinsic 
rest spectrum, so that in principle the counting could start 
at spectral type B3. 
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Table 4. Counts of OB main sequence stars with and with- 
out emission in the SMC per spectral type with OGLE-III. 



B2 B3 B4 B5 B6 B7 

Spectral types 



Fig. 4. Spectral-types distribution of OB stars in the SMC 
counted with OGLE-III. 



As a consequence of the many effects determining the 
apparent spectral type of Oe/Be stars, the mass-spectral 
type mismatch can be high. Since the number of stars in- 
creases rapidly as the spectral type becomes cooler, we pre- 
fer to deal with a lower limit of counted candidates. We thus 
considered stars only from the apparent spectral type B2 
and hotter. Once the number of all 0/B stars, i.e. with 
and without emission lines, that enter the chosen range of 
spectral types in the SMC, we determine the number of 
O e/Be stars using t he fre quencies of Be and Oe stars given 
bv lMartavan et al] (|201Cl[ ). which are based on a slit less sur- 
vey of emission-line stars in the SMC. 

The counts of all 0/B stars entering the desired spectral 
type ranges were perform ed using the new OGLE-III cat- 
alog (|Udalski et al.ll2008D . which covers 14 square degrees, 
so that both the bar and the external fields of the SMC arc 
covered. We note that OGLE-II data cover only 2.4 square 
degrees. Its use would require extrapolating the number of 
OB stars for regions off the bar to infer the total number 
of 0/B stars in the SMC with the required spectral types. 

The stars were selected using the dereddened color di- 
agram [My, {V — I)o]- To do the correction for the ISM 
extinction, we calculated the mean E[B-V] color excess 
based on the individual valu es from the SMC clusters (see 
iPietrzvnski &: Udalskil [19991) . which amounts t o 0.084 mag 
and is similar to the E[B-V]=0.08 mag given bv lKeller et al.l 
(|l999t) . Then the absolute My and the dereddened {V-I)o 
were obtained f o llowin g the same process as descri bed in 
iMartavan et al. (120101). Using the calibrations by iLand 
(|l992fl and lSaear fc CannonI (|l993), the E[V-I] is obtained 
by the formula: 



E[V - /] = 1.44 X E[B -V]+ 0.0175 x E[B - Vf 



(4) 



where the second term is negligible with E[B-V]=0.084. 
And finally the (V-I)o 



{V - /)o = {V-I)- 1.44 X E[B - V] 



(5) 



The number of OB stars in the SMC with and without 
emission lines so determined per subspectral type from B3 
to 08 is shown in Fig. 21 and detailed in Table H] for types 
between B3 and 08. The distribution in Fig. |4] shows that 
O and B-type stars basically follow the IMF except that 
the OGLE-III survey is not complete towards the late B 



08 


09 


BO 


Bl 


B2 


B3 


566 


671 


3482 


6269 


14406 


14332 



stars, which explains the sudden decrease in counts at B9. 
We are finally left with the following reference number of 
objects to predict the rate of LGBRs events: 



A^tot 

Mot 

Ntot 
Ntot 



4719 stars hotter than BO 
10988 stars hotter than Bl 
25395 stars hotter than B2 
39726 stars hotter than B3. 



(6) 



Since the FWHM of the OGLE-III PSF function is bet- 
ter than 1.34", it implies that in the SMC the counted 
early-type stars must be separated by at least 0.33 pc for 
their images not to merge. According to a recent study by 
ISabbi etaD (|2009l ). the highest early- type star density in 
the SMC is 8.5 star/pc^, which implies mean stellar sepa- 
rations of 0.34 pc. This ensures that the counting in Table S] 
is complete. 

In the next section we estimate the rates of LGRBs 
events and the number of LGRBs in the local universe. 
Their purpose is to get the range of uncertainties that can 
affect the estimate of rates depending on those inherent to 
stellar counts and the angles of beams. 

4.2. Prediction of the LGRB rates 

To predict the number of LGRBs seen from Earth, wc con- 
sider only specific beaming angles. In a first step, we de- 
termine the number of Oe/Be stars per spectral- type. We 
use for that the Be/ fB-|-Be) and Oe/ (0-l-O e) ratios in the 
SMC determined bv lMartavan et al.l (|2010t) . The numbers 
of stars being sought are given in column 2 of Table [51 
where the ranges stand for the uncertainties on the Be/B 
fractions. The yearly frequency of LGRB events is then ob- 
tained by dividing the counts per mass category by the 
respective total lifetime. The stellar ages used w ere interpo- 
lated among those obtained for this purpose bv lYoon et al.l 
(|2006D for stars between 13 and 30Mq with SMC metalhc- 
ity. Assuming that the LGRBs only refer to single fast- 
rotating staJSj_wB_j^move the b inaries from the counts. 
As from iPorter fc RiviniusI (|2003[ ) the binary rate is about 
30% among Oe/Be stars, the obtained above rates must 
be multiplied by 0.7; however, the exact ratio of binaries 
in B-type stars is far from being well known. Among the 
hottest ones and in particular f or O stars, the frequenc y of 
binary systems can reach 75% (|Sana et al.ll2008L I2009D . It 
is worth noting that only those binaries should be removed 
where the interaction between the components is signifi- 
cant, which is not always the case in Be-b i nary s ystems. 
Finally, we stress that the IMartavan et al.l (|2010f ) survey 
is a single-epoch survey. Because of t he variab l e char acter 
of the Be pheiiqmeno n, according to iFabreeatl (|2003[ ) and 
iMcSwain et al.l (|2008[ ) one third of the Be stars can be miss- 



ing. As a consequence, the final range of rates of LGRBs is 
calculated by multiplying the above rates by 0.7/(1-1/3). 

The LGRBs rates obtained for the SMC are called here 
the "LGRB base-rates" . They represent the intrinsic rates 
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that we generalize to all irregular galaxies of Magellanic 
type, which are metal poor and can host LGRBs progeni- 
tors. The LGRB base-rates arc thus given by 



R 



LGRB 



0.7 



1 - 1/3 



E 

Sp 



N, 



Oo/B( 



■{Sp) 



tiSp) 



LGRBs/yr(7) 



In this relation, A^oc/Bc(Sp) is the number of stars with 
spectral type Sp displaying the Be phenomenon, and t{Sp) 
is the total lifetime of stars with spectral type Sp in the 
early MS evolutionary phases. The constant factor in re- 
lation ([7]) is of the order 1, which indicates a fortuitous 
compensation between the fraction of rejected stars for bi- 
narity and the completeness factor for missing Oc/Bc stars. 
The -Rlgrb i'9'tes calculated with ([7]) are given in column 
3 of Tableil 

Since the LGRB phenomenon is strongly coUimated, to 
obtain the potential rates that can be observed, in a sec- 
ond step we have to multiply the base-rates by the the 
probability that beaming angle sweep t he Earth. The beam 



angles can ran ge from 0.5° to 15° (iHirschi et al I 



Frver et all 120071: iLamb et all 120051: IWatson et a" 



2005 



200a 



Zeh et aiT2006|) . The probability that a beam of angle 



be seen at Earth is: 



p{e) 



Jo 



sine de 



(8) 



Finally, to calculate the predicted yearly rate of LGRB 
events per average galaxy in a given volume of the Universe, 
we have to take the fraction of galaxies into account that 
have the required properties in this space to host LGRB 
progenitors. Let us then assume that the SMC properties 
regarding the frequency of fast rotators with low metallicity 
is partaken by the irregular Magellanic-type galaxies (Im- 
type), which are metal poor. We note this fraction as /im ~ 
Nim/NQ, where A^im is the number of Im-type galaxies in 
a given volume of the Universe, and A^g the number of all 
ga laxies in the same volume o f spac e. According to statistics 
by iRocca-Volmerange et all (|2007D . it is /ii„ = 0.17, which 
is a very lower limit and valid up to at least redshift z ~ 0.5. 

We present the results fo r the angles 9 = : 5 °, 10° , 
and 15°, which are used by iPodsiadlowski et al. (120041) 
Howev er, the angle 9 = : 5° according to I Watson et al.l 



^qM) is quite frequent. The angle 10° is roughly in the 
middle of the observed interval of opening angles and the 
angle 15° was chosen because is nearly the most opened one. 
although the least probable, but quite frequently evoked 
and used in the literature. It then happens that p = 0.0038 
for an angle of 5°, p = 0.015 for 10°, and p = 0.0341 for 15°. 
As there are 2 cones in each LGRB event, the probability 
([5]) has to be multiplied by a factor 2. Thus, the predicted 
yearly rates of LGRB events per average galaxy produced 
in the local Universe, ^^lq^b, given in columns 4 to 6 of 
Table [5] were obtained from 



<GtB(^) = ^:LGRBX2p(0)x/i„ 



(9) 



The observed yearly rate of LGRB events per aver- 
age galaxy that we use as a reference in this work to 
which we compare the predicted ones from OeBe stars 
is obtained from the total rate of GRBs, of which 2/3 
are considered of long-duration. From the BATSE mon- 
itoring it is found that i?£GRB ^ (0-2 - 3) x 10"'^ 



LGRBs/year/galaxy, where the extreme values concern 
the lo cal Un iverse or the whole Hubble volume, rcspcc- 
tivelv ll Zhang fc Meszaros 2004: PodsiadlQwski et al. 200l 
iFrver et al.E007t ). According to the (r2/»n: M /Afg ) -area of 
LGRBs progenitors indicated by Yoon et al.l (120061) and th e 
mass-calibration bv lHuang fc GiesT (|2006D or iLmS (IT991) . 
but neglecting spectral type changes from rotational and 
circumstellar effects discussed in Sect. SI we should only 
consider counts of stars hotter than spectral types BO-Bl. 
Taking the beaming angles from 5° to 15° into account, we 
note that the estimates made with populations of Be/Oe 
stars from BOe to 08e overlap the observed range of LGRBs 
rates adopted above. The hottest populations (09e to OSe) 
alone are not able to reproduce the rates. 

Another way of comparing the predicted rates with the 
observed ones is to determine the number of LGRBs ob- 
served in the local Universe up to a given redshift z. The 
number of LGRB events that have taken place in the local 
Universe and have been seen from Earth during the last Y 
years, is given by: 



-'^LGRB 



^ctB^X^GXy, 



(10) 



where Nq is the number of galaxies of all type counted 
up to a redshift z. We chose z = 0.2, since the number of 
galaxies in this volume see ms to be about well de termined. 
From the 2MASS catalog (jSkrutskie et all l2006fFI . we ob- 
tain Ng{z < 0.2) = 1140 931 galaxies. Adopting Y = 11 
years (see below), the predicted number of LGRB events 
seen from Earth are given in Table [6] 

Using the GRBox from the University of California at 
BerkelejO for years between 1998 and 2008, there are 8 
LGRBs (080108, 060614, 060505, 060218, 051109B, 031203, 
030329A, 980425) and 3 SGRBs (061201, 050709, 000607) 
with a redshift lower than 0.2. This proportion of 73% of 
LGRB and 27% of SG RB agrees with the propo rtion of 
75-25% as reported by IZhang fc Meszaros! (2Q04\ for the 
L GRBs-SGRB s. 

IZeh et all (|2006D and IWatson et cH (|2006[) have stud- 
ied the distribution of opening angles, which is reproduced 
in Fig. O If the LGRB event rate was estimated consider- 
ing this distribution of beam angles, the '2p{9) probability 
factor in relation ^ had to be replaced by the integrated 
probability P given by 



16° 



P = 2 X / p{9)Q{9)d9 = 0.015 



(11) 



where p{9) = 1 — cos 6* is given by ([8|), and Q{9) is the 
distribution shown in Fig. [SJ From ([TT]) we deduce that 
P defines an average opening angle of 7°, roughly. We can 
then interpolate the corresponding LGRB rate or number of 
events in Tables 5 and 6, respectively, which gives Nl'^'^'^ 



(2 - 5) X 10-^ LGRBs/yr/galaxy or Nf""^ 



LGRB 



'LGRB 

^3-6 LGRB 

events in the past 11 years. On account of the uncertainties 
that may affect not only our knowledge of the opening angle 
of beams or their distribution, but also the estimated rates 
of observed LGRB events in the local and Hubble Universe, 
we may consider that the counts limited to the spectral 



see also http://www.haydenplanetarium.org/universe/ 
duguide/exgg_twomass.php ' 

4 

see 



http://lyra.berkeley.edu/grbox/grbox.php7starttime 
=670702&endtlme=0§12-il 
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Table 5. LGRBs rates estimated with the first approach. 



Mass category 


Number of Bc/Oc stars 


LGRBs base rates 


Proba 5° 




Proba 10° 




Proba 15° 


B2e to 08e 


497S-6110 


3.9-4.S xlO"^ 


2.5-3.1 xlO- 


-7 


1.0-1.2 xlO 


-7 


2.3-2.S xlO-6 


Ble to OSe 


2774-3244 


2.6-3.0 xlO"'' 


1.7-1.9 xlO' 


-7 


6.7-7.7 xlO 


-7 


1.5-1.7 xlQ-^ 


BOe to OSe 


14S3-1551 


1.7-l.S xlO-4 


1.1-1.2 xlO" 


7 


4.4-4.6 xlO 


-7 


0.9S-1.0 xlO-6 


09e to OSe 


257-294 


5.3-6.2 xlO-5 


3.4-4.0 xlO- 


-8 


1.4-1.6 xlO 


-7 


3.1-3.6 xlO-^ 


OSe 


llS-135 


2.4-2.S xlQ-^ 


1.6-l.S xlO- 


-8 


6.2-7.2 xlO 


-8 


1.4-1.6 xlO-^ 



Table 6. Number of LGRBs predicted in the local universe (z<0.2) in 11 years. 



Mass category 


Number for anglc=5° 


Number for angle=10° 


Number for angie=15° 


B2e to OSe 


3-4 


11-14 


25-31 


Ble to OSe 


2-2 


7-9 


16-19 


BOe to OSe 


1-1 


5-5 


11-11 


09e to OSe 


0-0 


2-2 


3-4 


OSe 


0-0 


1-1 


2-2 



type interval BO-Bl to 08 predict rates that overlap rather 
well with the observed ones. This is the first time that an 
estimate of the LGRBs rates is done on the basis of a well- 
defined population of stars, i.e., born fast rotators with low 
initial metallicity. 

Very recently, iGeorgy et al.l (|2009t ) has tried to repro- 
duce the observed rates of LGRBs with models of WR stars. 
Although they obtained higher values of LGRBs rates than 
what observations seem to indicate, these authors conclude 
that only a fraction of WC /WO stars among those that ini- 
tially were rapid rotators could produce gamma ray events. 
This fast-rotation argument may then focus some interest 
on stars with the characteristics discussed in the present 
work. 

The counts of Oe/Be candidates for LGRB progenitors, 
as they were done here, however, represents a mere lower 
limit. Stars with masses corresponding to rest BO-Bl spec- 
tral types that rotate at fi/f^c 0.7 following the require- 
ments shown in Fig. |3] then display apparent spectral types 
cooler by at least one subspectral type. A prudent lower 
limit set at spectral type B2, would nevertheless enhance 
our predictions by a rough factor 2, as seen from counts 
displayed in (jH]), or Tables [5] and [5] However, massive fast 
rotators without emission lines should also be considered 
potential LGRBs progenitors as we see in Sect. 15.21 so that 
the final predicted rates may widely surpass the observed 
ones. 



5. Other phenomena and possible biases in the 
LGRBs rate determination 

5.1. Binaries 

If the LGRBs could also occur in a binary system 
(jCantiello et al.ll2007l: Ivan den Heuvel fc Yoonll2007D under 
the same conditions as fast rotation and chemical composi- 
tion, the rates we estimated abo ve should be multip l ied by 
1/0.7. However, the results from lLanger fc NormanI (|2006[ ) 
suggest that the binary evolutionary path is unlikely to 
produce LGRBs when they a re associated with metal-poor 
galaxies (jKewlev et al.ll2007f ). 




Fig. 5. Frequency of beam angles presented by 
IWatson et aLr(|2006f ) and lZeh et all (|2006D . 



5.2. Fast rotators without emission lines 

The ZAMS rotational velocities of Oe/Bc stars in the MW 
form the tail of a true rotational velocity distribution ^{V) 
that starts a t a ve locity limit given by Vl = 30 x (A/ /A/©) 
(jZorec et al.l l2007() . The tail of the velocity distribution 
formed by the Oe/Be stars can be represented with a 
Gaussian function, whose dispersion is also mas s-dependent 
and is roughly given by cry = 28 x (M/Mq) (IZorec et al.l 
|2007| ). This implies that there are rotators that can only 
acquire the Be phenomenon at the very end of their MS 
evolutionary phase, which cannot be seen as Oe/Be stars 
today. On the other hand, there are also fast rotating stars 
in the ZAMS that have slightly lower rotational velocities 
than Vi^ required to become Oe/Be, which will never dis- 
play emission lines. However, it can easily be shown that 
they have a global specific angular momentum that is equal 
to or larger than given by Eq. [2l so that they certainly have 
enough angular momentum to be considered as candidates 
for LGRBs progenitors. 
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Let us also take into account those objects which in 
principle have the required angular momentum but were 
not included in our counting to derive the LGRBs rate. 
To this end we assume that they are rigid rotators. The 
equatorial velocity of all those having the specific angular 
momentum given by Eq.[2]is readily calculated to give Vj = 
100 X (M/Mq)°-28. The fraction of objects yet not included 
in the estimated LGRB rates is then represented by the 
correcting factor q calculated as 



1 - erf(a:L) 
1 - erf (xsbh) 



(12) 



for which, as noted above, we have taken the high velocity 
tail ^{V) given by a Gaussian function. In Eq. [T^lerf (cc) is 

the error function, xl = Vl/v^ctvj and xsbh = V]/\/2ay. 
The factor q then ranges from 1.2 to 1.5 for masses from 
17 to 30Mq. We can then multiply the LGRB rates found 
in Sect. H] by an average factor g ~ 1.35, which increases 
even more the predicted rate well above the observed one 
and makes the massive fast-rotating stars highly probable 
progenitors of LGRBs. 

Finally, from the equatorial velocity Vj required to at- 
tain the stellar specific angular momentum of the order 
given by Eq. [2l we derive the angular velocity ratios that 
range from il/ilc = 0.50 to 0.52 for stars in the E-catcgory 
in the ZAMS. This means that the LGRBs progenitor area: 
E mass-category of Fig. |3l can be entirely covered by the 
0/B-type fast rotators. Assuming that the distribution of 
the rotational velocities of the stars e ntering the LGRBs 
progenitor region in Fig. [3] foreseen by lYoon et all (|2006D 
is the same as used to calculated the fraction q in relation 
(fT^ . it follows that our E mass group represents 44% of 
the stellar population in the LGRBs region. 

5.3. Differential rotators 

Models of stars that start in the ZAMS as rigid rotators, 
transform their original internal uniform angular veloc- 
ity into a profile of differential rotatio n, where the stel- 
lar cor e rotates faster than the envelope iMevnet fc Maedeil 
()2000t) . in no more than 2-3% of their MS lifespan. The 
latest theoretical works of stellar evolution with rotation 
claim, however, th at mag i ietic fi elds through the Tayler- 
Spru i t instability (ISpruiti Il999l 120021: iMaeder fc Mevnetl 
120051: lYoon et al.ll200fih are able to maintain the rigid ro- 
tation. Nevertheless, a recent critical dis cussion of the efii - 
ciency of the Taylcr-Spruit instability bv lZahn et al.l (|2007[ ) 
shows that the action of magnetic fields can be much less ef- 
ficient than previously thought. In all works of stellar struc- 
ture, it is also assumed that convection imposes f2 — con- 
stant. 

A number of discussions show that convection likely re- 
distributes the specific angular momentum j = Q,r^ rather 
than il, which implies that in the convection zones the an- 
gular vel ocity should be O oc 7-"^, w here p is a positive 
quantity (jTavledll973l: lDeupreelll998[) . In the Sun, the dif- 
ferential rotation exists precisely in the convective zones! 
Since the stars do not start their life in the ZAMS, but 
have passed though a long story before this phase, their 
internal rotation could already be differential in the ZAMS 
and with a larger content of angular momentum than just 
the allowed by a profile of rigid rotation. For an average 
star of 2OM0, it is then a simple matter to calculate what 



must be in order of magnitude its ZAMS ratio of energies 
T = K/\W\, higher than some maximum Tcrit, rigid = 0.015 
allowed by the rigid rotation {K is the rotational kinetic 
energy; W is the gravitational potential energy) in order 
for the star to keep up the same appearance as a rigid ro- 
tator throughout its MS evolutionary span. For that, we 
simply ask that the star fulfill at least two requirements: a) 
at any Lagrangian mass-coordinate the star never reaches 
the fateful limit t~0.14, which can make it secularly un- 
stable and would tend to transform the delimited stellar 
volume into a three-axial Jacobian ellipsoid; b) during the 



whole MS phase it has on the surface Vc 



equator 



critical- 



Calculations show that these conditions can be fulfilled by 
TzAMS — 0.03, which is roughly two times Tcrit, rigid in the 
ZAMS. 

If by chance O /B stars started their MS life span as neat 
differential rotators with r > Tcrit, rigid, the internal mixing 
processes can be strong in a large number of stars where 
this phenomenon is unexpected a priori because of their 
apparently slow rotation. Then make their chemical com- 
position attain the required properties to become potential 
LGRBs progenitors. Consequently, the estimate of the SMC 
-Rlgrb rate we made in this paper might still be increased 
with stars of a wider variety of masses and apparent low 
surface velocities. 



5.4. Underluminous GRBs 

Recently, it has been discovered tha t underluminous GRBs 
(iFrver et al.ll2007l: iFolev et al.ll2008D could not be observed 
by the gamma ray satellites or not be classified as a GRB. 
According to these authors, this finding indicates that there 
may be a population of events that are less luminous but 
are possibly from 10 to 10^ times more frequent. If these 
events are also considered, the local rate of GRBs could 
then amount to ~ 10~^ — 10~^ LGRBs/year/galaxy. If cor- 
rect, this last rate will indicate that other masses or kinds 
of stars can also be concerned such as less-massive B stars, 
in particular Bn stars, these stars are objects with spectral 
types that are most frequently later than B7. Moreover, 
they are fast rotators that do not have emission lines. 



5.5. Other possible progenitors 



According to iHammer et al.l (j2006[ ). runaway WR stars 
could be another potential progenitor, but it is very dif- 
ficult to estimate their impact on the LGRBs rates. 

5.6. The "first stars" and the LGRBs 

Using the GRBox from 1998 to 2009, we computed the 
statistics of GRBs by type with intermediate (z>3) and 
high redshifts (z>5). In the intermediate-z category, most 
of the GRBs are found to be long-GRBs (29 LGRBs, 2 
SGRBs, 8 unknown). In the high-z category, the result is 
similar, 6 are LGRBs or type II, 2 are unknown, and are 
SGRBs. The very high redshift GRB08091 3 was reclassified 
as a type II GRB by IZhang et aP (|2009l ). And the most 
distant object known, GRB 090423, whic h has a redshift 
of 8.2 , is also an L GRB or ty pe II GRB ([Salvaterra et al.l 
20091: iTanvir et all [200 9: Zh ang et al.ll2009D . According to 
Salvaterra et al. " (r2009[ ). the properties of the very- high red- 
shift LGRBs are similar to those of low/intermediate red- 
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shift. According to lYoon et al.l (j2006|) models, the number 
of LGRBs should increase with the redshift and the de- 
crease in Z until the first stars. However, the number of 
LGRBs, which can originate in Pop III stars, is limited 
by this type of population. At high redshift it is expected 
that the stars have very low metallicity and are similar to 
the first stars. In this case, they can reach rotational ve- 
loci ties as high as 8 km s~^, o r even higher, a s show n 
bv IChiappini et al.l (120061) and lEkstrom et all (l2008aD . 



lEkstrdm et al.l (|2008a^ and lMevnet fc Maeded (|2007[ ) show 
that these stars can also retain enough angular momentum 
in their core to c ollapse into bl ack holes and produce an 
LGRB. Moreover. iHirschil (|200l1 ) show that the stars with 
extremely low metallicity can rotate very fast and produce 
LGRBs. All these theoretical results agree with the statis- 
tics reported above and indicate that the stellar rotation 
close to the breakup velocity, as in the Be-phenomenon 
at low metallicity, is a key ingredient for the production 
of LGRBs. As a consequence, the LGRBs at high redshift 
should find their origin in the fast-rotating first stars. 



dictions by more than 30%. On account of all possible un- 
certainties that may affect the calculated and the observed 
rates of LGRBs in the local Universe, we can consider that 
both overlap correctly. Young star formation regions in the 
SMC like NGC346 could then have hosted or will host one 
LGRB in the next Myrs. 
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6. Conclusion 

While the first models focused primarily on fast-rotating 
WR stars as probable progenitors LGRBs, the latest devel- 
opments indicate that fast-rotating, low-metallicity massive 
B- and 0-type stars could also be behind the LGRBs. Such 
massive B- and 0-type stars can actually be the Be and 
Oe stars of low metallicity found in the SMC. To test this 
assumption, we used fundamental parameters derived for 
Oe/Be stars in the SMC observed with VLT-FLAMES in- 
struments. We first compared the average Vs'mi determined 
for several mass subsamples of SMC Oe/Be stars with the 
theoretically predicted rotational velocities. We found that 
today these objects on average have a surface angular ve- 
locity ratio Q/Qc— 0.99. Then, we determined the ZAMS 
rotational velocities again for Be stars in the SMC that we 
had studied in previous works and determined the ZAMS 
rotational velocities for a new subsample composed by SMC 
Oe stars. The average angular velocity rate of all these ob- 
jects in the ZAMS is (17/17c) =0.95. The ZAMS n/n^ rates 
thus determined were compared with those of fast-rotating 
massive stars foreseen by model calculations of LGRBs pro- 
genitors having an initial metallicity similar to the average 
one in the SMC. This comparison indicates that massive 
Oe/Be stars in the SMC could certainly have the proper- 
ties required to be taken as LGRBs progenitors. A lower 
limit of 17% of galaxies in the near Universe are of irregu- 
lar Magellanic-type that, by definition, all have metallicities 
Z < 0.002 - 0.004. AU can then host LGRBs progenitors. 
From the discussion it appears that not only can massive 
Oe/Be stars be LGRBs progenitors, but a larger popula- 
tion of 0/B stars without emission, which are less massive 
and with apparent lower surface rotation velocities, can also 
partake of this quality. The total number of LGRBs esti- 
mated from the Oe/Be star populations based on spotted 
beaming angles and counts from BO-Bl to 08 stars is on 
average 

-Rlgrb^IO"^ to -10-6 LGRBs/yr/galaxy, which 
represent -/Vlgrb ~ 2 to 14 in 11 years. These predictions 
increase beyond the observed ones, if the stellar counting 
starts at spectral types B1-B2, in accordance with the ap- 
parent spectral types of massive fast rotators. Massive fast 
rotators without emission lines can still enhance these pre- 



References 

Ballercau, D., Chauvillc, J., & Zorec, J. 1995, A&AS, 111, 423 
Bouret, J.-C, Lanz, T., Hillier, D. J., et al. 2003, ApJ, 595, 1182 
Campana, S., Panagia, N., Lazzati, D., et al. 2008, ApJ, 683, L9 
Cantiello, M., Yoon, S.-C, Langer, N., & Livio, M. 2007, A&A, 465, 
L29 

Chandrasekhar, S. & Munch, G. 1950, ApJ, 111, 142 
Chauville, J., Zorec, J., Ballereau, D., et al. 2001, A&A, 378, 861 
Chiappini, C, Hirschi, R., Meynet, C, et al. 2006, A&A, 449, L27 
Dcupree, R. G. 1998, ApJ, 499, 340 

Domiciano de Souza, A., Kervella, P., Jankov, S., et al. 2003, A&A, 
407, L47 

Ekstrom, S., Meynet, G., Chiappini, C, Hirschi, R., & Macder, A. 

2008a, A&A, 489, 685 
Ekstrom, S., Meynet, G., Macder, A., & Barblan, F. 2008b, A&A, 

478, 467 

Fabregat, J. 2003, in Astronomical Society of the Pacific Conference 
Series, Vol. 292, Interplay of Periodic, Cyclic and Stochastic 
Variability in Selected Areas of the H-R Diagram, ed. C. Sterken, 
65-+ 

Foley, S., McGlynn, S., Hanlon, L., McBreen, S., & McBreen, B. 2008, 
A&A, 484, 143 

Fremat, Y., Neiner, C, Hubert, A.-M., et al. 2006, A&A, 451, 1053 
Fremat, Y., Zorec, J., Hubert, A.-M., & Floquet, M. 2005, A&A, 440, 
305 

Fryer, C. L., Mazzali, P. A., Prochaska, J., et al. 2007, PASP, 119, 
1211 

Georgy, C, Meynet, G., Walder, R., Folini, D., & Macder, A. 2009, 
A&A, 502, 611 

Gouliermis, D. A., Chu, Y.-H., Hcnning, T., et al. 2008, ApJ, 688, 
1050 

Hammer, F., Flores, H., Schaerer, D., et al. 2006, A&A, 454, 103 
Hirschi, R. 2007, A&A, 461, 571 

Hirschi, R., Meynet, G., & Maeder, A. 2005, A&A, 443, 581 
Huang, W. & Gies, D. R. 2006, ApJ, 648, 580 

Hunter, I., Lennon, D. J., Dufton, P. L., et al. 2008, A&A, 479, 541 
Iwamoto, K., Mazzah, P. A., Nomoto, K., et al. 1998, Nature, 395, 
672 

Iwamoto, K., Nakamura, T., Nomoto, K., et al. 2000, ApJ, 534, 660 
Kawaler, S. D. 1987, PASP, 99, 1322 

Keller, S. C, Wood, P. R., & Bessell, M. S. 1999, A&AS, 134, 489 
Kewley, L. J., Brown, W. R., Geller, M. J., Kenyon, S. J., & Kurtz, 

M. J. 2007, AJ, 133, 882 
Lamb, D. Q., Donaghy, T. Q., & Graziani, C. 2005, ApJ, 620, 355 
Lang, K. R. 1992, Astrophysical data, Planets and Stars, New York, 

Springer Verlag Eds 
Langer, N. & Norman, C. A. 2006, ApJ, 638, L63 
MacFadyen, A. I. & Woosley, S. E. 1999, ApJ, 524, 262 
Maeder, A., Grebel, E. K., & Mermilliod, J.-C. 1999, A&A, 346, 459 
Maeder, A. & Meynet, G. 2001, A&A, 373, 555 
Maeder, A. & Meynet, G. 2005, A&A, 440, 1041 



12 



Christophe Martayan et al.: SMC Be/Oe stars as LGRBs progenitors. 



Martayan, C, Baade, D., & Fabregat, J. 2010, A&A, 509, All 
Martayan, C, Fremat, Y., Hubert, A.-M., et al. 2007, A&A, 462, 683 
Martins, F., Hillier, D. J., Bouret, J. C, et al. 2009, A&A, 495, 257 
McSwain, M. V., Huang, W., Gies, D. R., Grundstrom, E. D., & 

Townsend, R. H. D. 2008, ApJ, 672, 590 
Meilland, A., Stee, P., Vannier, M., et al. 2007, A&A, 464, 59 
Meynet, G. & Maeder, A. 2000, A&A, 361, 101 
Meynet, G. & Maeder, A. 2002, A&A, 390, 561 
Meynet, G. & Maeder, A. 2007, A&A, 464, Lll 

Modjaz, M., Kewley, L., Kirshner, R. P., et al. 2008, AJ, 135, 1136 
Moujtahid, A., Zorcc, J., & Hubert, A. M. 1999, A&A, 349, 151 
Moujtahid, A., Zorec, J., Hubert, A. M., Garcia, A., & Burki, G. 1998, 

A&AS, 129, 289 
Nota, A., Sirianni, M., Sabbi, E., et al. 2006, ApJ, 640, L29 
Pasquini, L., Avila, G., Blecha, A., et al. 2002, The Messenger, 110, 1 
Pietrzynski, G. & Udalski, A. 1999, Acta Astronomica, 49, 157 
Podsiadlowski, P., Mazzali, P. A., Nomoto, K., Lazzati, D., & 

Cappellaro, E. 2004, ApJ, 607, L17 
Porter, J. M. & Rivinius, T. 2003, PASP, 115, 1153 
Rocca-Volmerange, B., de Lapparent, V., Seymour, N., & Fioc, M. 

2007, A&A, 475, 801 
Sabbi, E., Gallagher, J. S., Tosi, M., et al. 2009, ApJ, 703, 721 
Sagar, R. & Cannon, R. D. 1997, A&AS, 122, 9 

Salvaterra, R., Vallc, M. D., Campana, S., ct al. 2009, Nature, 461, 
1258 

Sana, H., Gosset, E., & Evans, C. J. 2009, MNRAS, 1653 
Sana, H., Gosset, E., Naze, Y., Rauw, G., & Linder, N. 2008, MNRAS, 
386, 447 

Skrutskio, M. F., Cutri, R. M., Stiening, R., ct al. 2006, AJ, 131, 1163 
Spruit, H. C. 1999, A&A, 349, 189 
Spruit, H. C. 2002, A&A, 381, 923 

Tanvir, N. R., Fox, D. B., Levan, A. J., et al. 2009, Nature, 461, 1254 
Tayler, R. J. 1973, MNRAS, 165, 39 

Thone, C. C., Fynbo, J. P. U., Ostlin, G., et al. 2008, ApJ, 676, 1151 
Townsend, R. H. D., Owocki, S. P., & Howarth, I. D. 2004, MNRAS, 
350, 189 

Tutukov, A. V. & Fedorova, A. V. 2007, Astronomy Reports, 51, 847 
Udalski, A., Soszynski, I., Szymaiiski, M. K., et al. 2008, Acta 

Astronomica, 58, 329 
van den Heuvel, E. P. J. & Yoon, S.-C. 2007, Ap&SS, 311, 177 
van Marie, A. J., Langer, N., Yoon, S.-C., & Garcia-Segura, G. 2008, 

A&A, 478, 769 

Vinicius, M. M. F., Zorec, J., Leister, N. V., & Levenhagen, R. S. 

2006, A&A, 446, 643 

Watson, D., Hjorth, J., Jakobsson, P., et al. 2006, A&A, 454, L123 
Wisniewski, J. P. & Bjorkman, K. S. 2006, ApJ, 652, 458 
Wisniewski, J. P., Bjorkman, K. S., Magalhaes, A. M., et al. 2007, 

ApJ, 671, 2040 
Woosley, S. E. 1993, ApJ, 405, 273 
Woosley, S. E. & Bloom, J. S. 2006, ARA&A, 44, 507 
Woosley, S. E. & Heger, A. 2006, ApJ, 637, 914 
Yoon, S.-C, Langer, N., & Norman, C. 2006, A&A, 460, 199 
Zahn, J.-P., Brun, A. S., & Mathis, S. 2007, A&A, 474, 145 
Zeh, A., Klose, S., & Kann, D. A. 2006, ApJ, 637, 889 
Zhang, B. & Meszaros, P. 2004, International Journal of Modern 

Physics A, 19, 2385 
Zhang, B., Zhang, B., VirgiU, F. J., et al. 2009, ApJ, 703, 1696 
Zorec, J. 1986, Structure et rotation differentielle dans le etoiles B 

avcc ct sans emission (Paris; Universite VII, 1986) 
Zorcc, J., Ballereau, D., Chauville, J., & Garcia, A. 1992, Academic 

des Science Paris Comptcs Rcndus Scric B Sciences Physiques, 314, 

1429 

Zorec, J., Fremat, Y., & Cidalc, L. 2005, A&A, 441, 235 

Zorec, J., Fremat, Y., Martayan, C, Cidale, L. S., & Torres, A. F. 

2007, in Astronomical Society of the Pacific Conference Series, Vol. 
361, Active OB-Stars: Laboratories for Stellare and Circumstellar 
Physics, ed. A. T. Okazaki, S. P. Owocki, & S. Stefl, 539—1- 

Zorec, J., Mochkovitch, R., & Divan, L. 1988, Academic des Science 
Paris Comptcs Rendus Scric B Sciences Physiques, 306, 1265 



